Extrema of discrete Wigner functions and applications 
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We study the class of discrete Wigner functions proposed by Gibbons et al. [Phys. Rev. A 
70, 062101 (2004)] to describe quantum states using a discrete phase-space based on finite fields. 
We find the extrema of such functions for small Hilbert space dimensions, and present a quantum 
information application: a construction of quantum random access codes. These are constructed 
using the complete set of phase-space point operators to find encoding states and to obtain the 
codes' average success rates for Hilbert space dimensions 2, 3, 4, 5, 7 and 8. 



PACS numbers: 03.67.-a, 03.67.Ac 
I. INTRODUCTION 



The Wigner function W(q,p) was introduced by 
Wigner in 1932 [l| as a way to represent quantum states 
of one or more particles in phase-space. It is a quasi- 
probability distribution, which means it retains some of 
the properties of a true probability distribution, while 
having some surprising properties due to quantum ef- 
fects. For example, it can be negative in some regions in 
phase-space. 

There were many proposals of analogues of W(q,p) 
to represent quantum systems with discrete degrees of 
freedom such as spins (see 0] and references therein for 
a review). These discrete Wigner functions have been 
applied to visualize quantum states and operations in 
the context of quantum information and computation Q 
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In this paper we study properties and applications of 
the class of discrete Wigner functions defined by Gibbons 
et al. 0, which take values on a discrete phase-space 
built with finite fields. We start in sections HT1 and ITTTl by 
reviewing the definition of this class of functions. In sec- 
tion IIVI we calculate the spectra of the phase-space point 
operators used to define the discrete Wigner functions. In 
section|V]we describe how to use the spectra calculated to 
find the extremal values for the discrete Wigner function 
for some small Hilbert space dimensions. We also de- 
scribe how the phase-space point operators can be used 
in a quantum information application known as quantum 
random access codes. In section I VI I we introduce these 
codes with a simple example, and present a quantum 
random access code construction based on states that 
maximize the discrete Wigner function. 



II. DEFINING A CLASS OF DISCRETE 
WIGNER FUNCTIONS 

The discrete phase-space is a d x d grid in which we 
identify some particular sets of d points called lines. Par- 
allel lines are lines sharing no points in common. Follow- 
ing Gibbons et al. 0, a partition of the d 2 phase-space 
points into d parallel lines of d points each will be called a 
striation. The definition of lines and striations is done in 
such a way as to ensure, in this discrete geometry, some 
geometrical properties akin to the properties of lines in 
usual geometry: 

i given any two points, exactly one line contains both 

points; 

ii given a point a and a line A not containing a, there is 

exactly one line parallel to A that contains a; 

iii two non-parallel lines intersect at exactly one point. 

Gibbons et al. described how to define d(d + 1) lines, 
partitioned into d + 1 striations of d lines each, satisfying 
the requirements above. The construction is based on 
considering the discrete phase-space as a 2-dimensional 
vector space labelled by finite fields (for details, see (H). 
In [fj Wootters discusses different geometrical problems 
associated with this construction. 

To define a discrete Wigner function, we need to as- 
sociate a projector onto a quantum state to each line in 
discrete phase-space. These will be projectors onto a set 
of d + 1 mutually unbiased bases (MUB). Consider two 
different orthonormal bases B\ and B 2 : 

Bi = , \vi, 2 ) , ■-, \vi.d)}, |(vmKj}| 2 = Si,j, (1) 

B 2 = {\v2,l) , \V2,2) , \v2,d)}, \(v2,i\v2j)\ 2 = (2) 

These two bases B\ and B 2 are mutually unbiased if 

\(vi,j\vk,l)\ 2 = \ if i k. (3) 
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Wootters and Fields showed that one can define (d+1) 
such mutually unbiased bases for power-of-prime dimen- 
sion d [7j. Mutually unbiased bases have been studied 
because of their use in a number of quantum informa- 
tion applications, for example quantum cryptography [8|], 
quantum state and process tomography [9| and in the 
construction of quantum t-designs [lOj , used to estimate 
averages of functions over quantum states. 

To define a discrete Wigner function (DWF) we pick a 
one-to-one map between the lines in discrete phase-space 
and the projectors onto a complete set of MUB in the 
following way: 

• each basis set Bi is associated with one striation 

Si ; 

• each basis vector projector Qi j = (vij\ is as- 
sociated with a line A^j (the jth line of the ith 
striation). 

These maps define uniquely the values of the DWF Wp 
for all points (3 if we impose the following constraints: 

Tr(Q itj p)= ]T W a , (4) 

where p is the system's density matrix, and the sum is 
over phase-space points a in the line A associated with 
projector Qij. These requirements amount to demand- 
ing that the sum of the Wigner function over each line 
must be equal to the probability of projecting onto the 
basis vector associated with that line. 

Note that there are multiple ways of making these asso- 
ciations. In general, this will lead to different definitions 
of the DWF using the same fixed set of MUB. The proce- 
dure outlined above then leads not to a single definition 
of W, but to a class of Wigner functions instead. 

We now define the phase-space point operator Ap as- 
sociated with phase-space point f3: 

A p = Y,Qx-I, (5) 

AD/3 

where the sum is over projectors Q\ associated to lines 
A containing point f3, and / represents identity. The op- 
erators A/3 appear naturally when we invert equations 
(J4j) to write an expression for Wp in terms of the MUB 
projectors: 

W = - d Tr{pAp), (6) 

We see that the expectation value of Ap (multiplied by 
1/d) gives the value of the DWF at phase-space point f}. 

The A operators form a complete basis in the space 
of d x d matrices. It can be shown that the DWF value 
at point a is simply the expansion coefficient of p corre- 
sponding to the A a operator: 

p = Y / W a A a , (7) 

a 



The multiple ways of associating projectors with lines in 
phase-space result in multiple definitions for A a . While 
a single definition of W requires d 2 operators A a (one 
for each phase-space point), the full set of A a one can 
define with the same fixed complete set of MUB has d d+1 
elements. In section lVIl we will make use of the full set of 
phase-space point operators to obtain a construction for 
a quantum information application known as quantum 
random access codes. 



A. Negativity and non-classicality 

Cormick et al. [ll| have characterized the set of states 
which have non-negative Wigner functions. These states 
turn out to have some interesting properties, which we 
will review here, as they motivate the work reported in 
the remainder of this paper. 

For a (i-dimensional quantum system, one can find 
complete sets of (d+1) mutually unbiased bases using 
the finite-field construction introduced in (2j. This con- 
struction is only valid for power-of-prime d, since this is 
the necessary condition for a finite field to exist. One can 
then define the set Cd of e?-dimensional states which have 
non-negative DWF in all definitions, that is, whose ex- 
pectation values for all phase-space point-operators are 
non-negative. 

In UJ it was shown that the only pure states in Cd 
are the MUB projectors, which can always be chosen to 
be stabilizer states. The stabilizer formalism [l2| pro- 
vides a way to represent pure states in Cd using a num- 
ber of bits which is polynomial in the number of qubits. 
Since general pure states require a description which is 
exponential-size, the set Cd is classical in the sense of 
having a short description. 

For systems of prime dimensions, the two notions of 
classicality exactly coincide: the only pure states with 
nonnegative DWF are exactly the stabilizer states. In 
this context, negativity of any DWF (as witnessed by 
negativity of one of the A a operators) indicates non- 
classicality in the sense of the absence of an efficient de- 
scription using the stabilizer formalism. These results 
motivated us to investigate the extrema of the discrete 
Wigner functions. 

Before proceeding, we need to review some construc- 
tions of complete sets of MUB for rf-dimensional systems, 
as these are necessary to define DWF using eqs. ([5]) and 
©■ 



III. COMPLETE SETS OF MUTUALLY 
UNBIASED BASES 

For prime dimension d there is a canonical construc- 
tion of a complete set of (d + 1) MUB's first proposed 
by Ivanovic To review this construction, let |tv,k)j 
denote the j-th component of the A:-th vector in the r-th 
basis, r = 0,1,..., d. The vectors in the complete set of 
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TABLE I: Set of 5 MUB for two qubits 



TABLE II: Set of 9 MUB for three qubits 



1 


XX 


XI 


IX 


2 


zz 


Zl 


1Z 


3 


YY 


Yl 


1Y 


4 


XY 


YZ 


ZX 


5 


XZ 


YX 


ZY 



d+l MUB are then: 

Kfc)j 



\V r ,k) 



5jk 
1 



_L e w+i*) 



Vd 



1-3 A) 



|«(d-l),fc)j 
\Vd,k). 



1 

1 

— p« 



^((rf-l)j 2 +3fe) 



(8) 
(9) 

(10) 
(11) 

(12) 
(13) 

(14) 



When the Hilbert space dimension is a prime power, 
there are different constructions of complete sets of MUB. 
Let us now review a simple construction of a complete set 
of MUB for n qubits (Hilbert space dimension d = 2 n ) 
consisting solely of stabilizer states [14|. We start by 
considering the 4™ Pauli operators for n qubits, which 
are the tensor products of single-qubit Pauli operators 
X, Y, Z and identity. From this set, remove the iden- 
tity operator. The remaining 4™ — 1 Pauli operators can 
be partitioned into 2™ + 1 sets, each containing 2™ — 1 
mutually commuting Pauli operators. It was proven in 
fl4| that the common eigenstates of the operators in each 
such set form a basis, and moreover that the 2™ + 1 bases 
thus defined are mutually unbiased. 

We will now provide two examples of this construction, 
which will be useful to us later on. The first example is 
a set of 5 MUB for two qubits, each basis being formed 
by the common eigeinstates of each row of operators in 
Table HI In the table, operator XY for example stands 
for the tensor product of X on the first qubit by Y on 
the second. 

The second example is a set of 9 MUB for three qubits, 
comprising the common eigenstates of the operators in 
each row of Table [TTJ 



IV. SPECTRA OF PHASE-SPACE POINT 
OPERATORS 

As we have seen in the previous sections, the DWF is 
defined by eq. © using the phase-space point operators 



1 


XXX 


XXI 


XIX 


Xll 


1XX 


XXX 


XXX 


2 


XXY 


XYX 


YXX 


YYY 


ZZ1 


ZXZ 


XZZ 


3 


XXZ 


XYY 


YZ1 


YXX 


ZXY 


ZYZ 


XZX 


4 


XYZ 


XZX 


YXl 


Y1Y 


ZYX 


ZZZ 


XXY 


5 


XY1 


X1Z 


YXY 


YZX 


ZXX 


ZZY 


XYZ 


6 


XZY 


X1Y 


YZZ 


Y1Z 


ZZX 


ZXX 


XZX 


7 


XZZ 


XZ1 


YYZ 


XY1 


ZXZ 


ZXX 


XXZ 


8 


YXZ 


YYX 


YZY 


Yll 


XXZ 


XYX 


XZY 


9 


ZYY 


ZY1 


Z1Y 


ZU 


XYY 


XYX 


XXY 



A. In this section we calculate the spectra of A for the 
constructions of complete sets of MUB we reviewed in 
section Mil The spectra we tabulate in Table Mil agree 
with the spectra calculated independently by Appleby et 
al. [151 ] for Hilbert space dimensions d = 3, 4, 5 and 7. In 
Table UTT1 we also report the number of phase-space point 
operators with each spectrum, an information which will 
be necessary for the quantum information application de- 
scribed in section IVTl 

We have calculated the full spectra of all d d+1 phase- 
space point operators for d = 2,3,4,5,7 and 8, but the 
latter two cases have too many different spectra for us to 
reproduce here. We would like to point out only one piece 
of information about these cases: the extremal eigenval- 
ues found over all the A a . For d = 7 the largest eigen- 
value is X m ax — 2.4178 and the smallest is \ m in = — 1, 
whereas for d = 8 the largest is X ma x — 2.5490 and the 
smallest is X min = —0.9979. 

For the case of two qubits (d = 4) we calculated the 
spectra for all six different stabilizer MUB constructions 
of the kind described in [l4| , and the spectra found are 
identical to those of the MUB set in TablelH We have also 
established that for three qubits (d = 8) , one can find 960 
different stabilizer constructions of this kind. Testing a 
few of those we found exactly the same spectra as for the 
example in Table ITT1 Based on this, we conjecture that 
the phase-space point operator spectra is independent of 
which MUB set construction one uses. 



V. EXTREMA OF DISCRETE WIGNER 
FUNCTIONS 

Unlike probability distributions, we have seen that 
the discrete Wigner function can assume negative val- 
ues. What are the extremal values that it can attain? 
Wootters conjectured that the minimal value min(W) 
that the discrete Wigner function could assume would 
be —1/d for odd-prime Hilbert space dimension d 
He also showed that min(W) = —0.183 for d = 2 



16]. 
17], 



and together with Sussman [18| found a maximum value 
of max(W) = 0.319 for some particular definitions of 
discrete Wigner functions for d = 8. In this section we 
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TABLE III: Spectra of phase-space point operators 



(1 


Number 


Spectrum 


2 


8 


r l , V3 l V3i 

^ 2 , 2 2 J 


3 


9 






{-1,1,1} 




72 






rn 1 i 1 V5\ 
l u , 2 ~ 2 > 2 2 J 




320 




{-0.50000, -0.50000, 0.13397, 1.86603} 


1 

4: 


320 




{-0.86603, -0.50000, 0.86603, 1.50000} 




384 




{-0.89680, -0.14204, 0.27877, 1.76007} 




1000 


{" 


0.70281, - 


-0.61803, -0.13294, 0.48666, 1.96712} 




2000 


{ 


-0.79859, 


-0.36221, 0.00000, 0.10661, 2.05419} 




2000 


{ 


-0.83607, 


-0.81000, 0.00000, 1.05469, 1.59139} 


5 


3000 


{" 


0.83726, - 


-0.58152, -0.09576, 0.62870, 1.88584} 




1000 


{" 


0.90039, - 


-0.64018, -0.14531, 1.06785, 1.61803} 




3000 


{ 


-0.90932, 


-0.48701, 0.00000, 0.46853, 1.92780} 




3000 




0.94658, - 


-0.51690, -0.18438, 0.93842, 1.70944} 




600 


{ 


-1.00000, 


-0.61803, 0.00000, 1.00000, 1.61803} 




25 


{ 


-1.00000, 


-1.00000, 1.00000, 1.00000, 1.00000} 



describe a general method for finding the extrema among 
all discrete Wigner functions definable with a fixed com- 
plete set of MUB, and use it to explicitly calculate the 
extrema for d = 2, 3, 4, 5, 7 and 8. 

Recall that the phase-space point operator A a associ- 
ated with the phase-space point a is given as the sum over 
MUB projectors associated with all phase-space lines 
that contain a: 

A a = ^2 Qx — I, (15) 

where the sum is over projectors associated with lines A 
containing point a, and / represents identity. 

Given a phase-space point operator A a , we want to 
find the minimum of its expectation value 

min (A a ) = min (tv l^> W J J ~ L ( 16 ) 

The minimum results when \ip) is the eigenvector asso- 
ciated with the smallest eigenvalue \ m in of A a (see the 
Appendix for a proof). We can evaluate the spectrum of 
A a to find its smallest eigenvalue A^ in . Then, using eq. 
((6]) , the most negative value for the discrete WF at point 
a will be given by 

min(W a ) = -min (A a ) = ^Knin- (17) 

To find the most negative value for the function, one 
needs to find all eigenvalues of all possible phase-space 
point operators. For a d-dimensional system there are 
d d+1 different phase-space point operators, only d 2 of 
which appear in any single definition of a DWF. Mini- 
mizing over a, we can use eq. (|17[) to obtain the smallest 
value that the function can attain. 



TABLE IV: Extremal values for DWF 



d 






2 


i(l + \/3) ~ 0.683 


i(l- v/3) ^ -0.183 


3 


+ ~ 0.539 


i 

3 


4 


0.4665 


-0.2242 


5 


0.411 


i 

5 


7 


0.3454 


1 
7 


8 


0.3186 


-0.1247 



The same reasoning can be used to obtain the maxima 
of the DWF, using the largest eigenvalue of any of the A a . 
Using the spectra tabulated in section llVl we obtained the 
extremal values of the DWF for small dimensions d, listed 
in Table HVl The results support Wootters' conjecture for 
odd-prime d. 



VI. AN APPLICATION: QUANTUM RANDOM 
ACCESS CODES 

In this section we review the quantum information 
protocol known as quantum random access codes, and 
present a code construction that relies on states max- 
imizing the discrete Wigner function. Let us start by 
recalling what these codes are using a simple example. 

Imagine a situation in which Alice encodes to classical 
bits into n bits (m > n), which she sends to Bob, who 
will need to know that value of a single bit (out of the 
to possible ones) with a probability of at least p. We 
may represent such an encoding/decoding scheme by the 
notation: to — > n. 

Prior to sending the n-bit message, however, Alice does 
not know which of the to bits Bob will need to read out. 
To maximise the least probability of success p, Alice and 
Bob need to agree on the use of a particular, efficient 
to — ► n encoding. 

We can consider the quantum generalization of this 
situation, in which Alice can send Bob n qubits of com- 
munication, instead of n bits. The idea behind these so- 
called quantum random access codes (QRACs) is very old 
by quantum information standards; it appeared in a pa- 
per written circa 1970 and published in 1983 by Stephen 
Wiesner 0. 

These codes were re-discovered in [2(|, where the ex- 
plicit comparison with classical codes was made. 



A. Example: 3 — * 1 QRAC with a qubit 

Let us illustrate the idea with a 3 — > 1 quantum ran- 
dom access code (QRAC) that encodes three bits into a 
single qubit. This QRAC was attributed to Isaac Chuang 
in ref. [2$. 
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Instead of concerning ourselves with the least decoding 
probability of success p, we will use as figure of merit 
the average probability of success p q . With three bits, 
Alice has 2 3 = 8 possible bit-strings &o&i^2- For each 
possibility she will prepare one particular state from the 
set depicted in Figure [T] These states lie on the vertices 
of a cube inscribed within the Bloch sphere, which is the 
representation of one-qubit pure states using spherical- 
coordinate angles 0, <fi: 



\i/j(e, </>)) = cos(0/2) |0) + exp(^) sin(0/2) |1) . (18) 



If Bob wants to read out bit bo he measures along the re- 
axis and associates a positive result with bo = 0. To read 
bits b\ [62] he measures along the y-axis [z-axis] and again 
associates a positive result with 61 = [62 = 0]. It is easy 
to see that Bob's average probability of success is given 
by p q = cos 2 (0/2) = 1/2 + v/(3)/6 ~ 0.79, where the 
angle is given in the caption to Figure [TJ The optimal 
classical 3 — ► 1 random access code succeeds only with 
average probability p c — 0.75, as can be checked easily 
through a search over all deterministic protocols. 

Note that the QRAC just presented uses decoding mea- 
surements which are projections onto the canonical set of 
MUB for a qubit, that is, the X, Y and Z bases. The en- 
coding states are found by optimizing the probability of 
correctly decoding each coding state. 




FIG. 1: Encoding states for the 3 — > 1 QRAC using a single 
qubit. Alice prepares one out of eight states on the vertices 
of a cube inscribed within the Bloch sphere, depending on 
her three-bit string. The angle 9 is such that cos 2 (9/2) = 
l/2+\/(3) /6 ~ 0.79, which is the probability of Bob correctly 
decoding a single bit out of the three. 



B. A QRAC construction from A a 

The full set of d d+1 phase-space point operators A a can 
be used to build a particularly symmetric set of quan- 
tum random access codes. The encoding states will be 
those maximizing each (A a ). Our results from section fVl 
showed that those are the largest-eigenvalue eigenstates 
of each A a . 

Our goal is to construct a QRAC that encodes (d+1) 
messages with d possible values each, using a single quan- 
tum d-level system that will be sent by Alice and mea- 
sured by Bob. As in the case with qubits, our strategy 
will be for Bob to decode by performing projective mea- 
surements onto one out of the (d+1) MUB that exist for 
this d-dimensional system (for power-of-prime d). Alice 
has to find d d+1 different encoding states, each of which 
will decode correctly with the highest possible probabil- 
ity. 

As we have shown in sectionfVl the pure quantum state 
\t/j a ) that maximizes (A a ) is the eigenvector with largest 
eigenvalue X^ ax of (A a ). Bob's decoding procedure in- 
volves measuring the encoding state \ip a ) onto one of the 
(d + 1) MUB's. His average probability of successfully 
decoding will be the average ^xj ^2\ Da (Q\)- To max- 
imize this decoding probability for encoding state \ip a )i 
we need to maximize 



max 



d+1 



max((A a ) + 1)(19) 



(\ a max + 1). (20) 



The average performance of this QRAC protocol can be 
found by averaging the probability of success of the d d+1 
encoding states, each corresponding to one phase-space 
point operator A a : 



1 



Pq = 



d d ^ 



!)■ 



(21) 



We see that the protocol's average success rate p q de- 
pends on the value of the sum of the largest eigenvalues 
of each of the d d+1 possible phase-space point operators 
A a . 

Using the MUB constructions described in section IIIII 
and the point operator spectra calculated in section HVl 
we were able to compute p q for a (d+1) — ► d QRAC using 
systems of dimension d = 2, 3, 4, 5, 7 and 8. The results 
are summarized in Table [V] The construction recovers 
the known success rate of the 3 — * 1 QRAC with a qubit. 

It is clear that this construction can be extended to 
higher-dimensional systems, provided a DWF can be de- 
fined for them. This will be the case for power-of-prime 
dimensions d, using for example the finite-field construc- 
tion given in 

It would be interesting to investigate how these pro- 
tocols fare against the optimal classical protocols. This 
would require evaluating the optimal probability of suc- 
cess for a (d+1) — ► d <i-level classical random access code, 
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TABLE V: QRAC success rate p q 



Dimension 


Pq 


d = 2 


|(3 + V3) ~ 0.789 


d= 3 


£ + j£ ~ 0.637 


d = A 


0.5424 


d = 5 


0.4700 


d= 7 


0.3720 


d = 8 


0.3372 



something that to our knowledge has not been done for 
d > 2. In [2(| some asymptotic results were obtained 
for large d, indicating that while there may be an advan- 
tage of quantum over classical for small dimensions d, 
this advantage practically disappears in the asymptotic 
regime. 

VII. CONCLUSION 

We have reviewed the definition of discrete Wigner 
functions given in used to describe quantum sys- 
tems in a discrete phase-space. We calculated the spec- 
tra of phase-space point operators for small dimensions, 
and used them to obtain the extrema of the discrete 
Wigner function for systems in Hilbert-space dimensions 
d = 2,3,4,5,7 and 8. We then described the protocol 
known as quantum random access codes (QRAC), and 
used the phase-space point operators to find encoding 
states and to obtain the efficiency of a QRAC construc- 
tion whose encoding states maximize the discrete Wigner 
function at each phase-space point. 



are obtained when \<fi) is an eigenstate corresponding to 
extremal eigenvalues of p. 

We start by showing that the state that maximizes this 
expression is the eigenstate \\ ma x) corresponding to the 
largest eigenvalue X ma x of p. It is easy to evaluate eq. 
for |0) = \X max ): 

(Xmax\p\ X m ax) X m ax (Xmax\Xmax) X 7n ax- (A'2) 

Let us now consider the expansion of a general state \4>) = 
Y^idi I*) m the basis that diagonalizes p. For \<j>), the 
expression we are trying to maximize takes the value: 

(<f>\p\<f>) =^A 4 |d l | 2 = A ma:c (l-A)+ A *l d *| 2 ' 

(A.3) 

where Xi are the eigenvalues of p. We have rewritten the 
sum to single out the term corresponding to X max , and 
defined A > such that (1 — A) = |rf mQ2; | 2 . Now, because 
all Xi < Xmax (by definition of X ma x), we have 

(4>\p\ (j)) = X m ax a-A)+ a *ki 2 < a 

max ■ 

(A.4) 

So we have proven that any set of coefficients different 
from those of \X m ax) leads to a smaller expectation value 
for p. and so |A maa; ) maximizes this expectation value. 
A similar argument can be made to prove that the state 
that minimizes (4>\p\ 4>) is the eigenstate corresponding to 
the smallest eigenvalue of p. 

The claim following eq. (|16p is justified by applying 
the result above to p = -^-^ ^2\ Da Q\- 



APPENDIX 

Let p be a density matrix, hence positive semi-definite, 
Hermitian and with unit trace. Let \tp) be a pure state. 
In this Appendix we prove that the extrema of 

m<t>) (a.i) 
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